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Abstract 



The color-evaporation model quantitatively describes all data on photo- 
production and hadroproduction of charmonium. Although the model is in 
part nonperturbative, the associated parameters can for instance be deter- 
mined from the charmonium photoproduction data. At this point its predic- 
tions for the prompt production of t/^'s at the Z pole are made with no free 
parameters. We show here that this approach successfully describes all data 
on the production of prompt -i/^'s in Z decays. 
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I. INTRODUCTION 



There has been a renewed interest in studying the mechanism by which charmonium 
is produced, triggered mostly by some puzzhng data from the Fermilab Tevatron. The 
standard perturbative QCD calculations using the color-singlet model failed to explain the 
data, occasionally by orders of magnitude ||T|. The data incited a complete review of the 
treatment of color in QCD and can, in fact, be explained by allowing perturbative color 
octet cc states to evolve into the asymptotic colorless charmonium states. This prescription 
is present in both the color-evaporation and in the color-octet models 0. 

At the Z resonance the majority of ■^/''s are produced via 6-hadron decays, with a branch- 
ing ratio B(Z ^ + X) of approximately 4.0 x 10^'^ 0. Recently the OPAL collaboration 
also presented evidence for the prompt production of il) in hadronic Z decays 

We study here the prompt charmonium production in Z decays in the framework of the 
color-evaporation model. The predictive power of the model stems from the fact that its free, 
nonperturbative parameters can be fixed by charmonium photoproduction data, which leads 
to very well defined predictions for the inclusive charmonium branching ratios. Although 
part of the mechanism is essentially nonperturbative, it can be described by a minimal set 
of free parameters which are few in comparison with the color-octet model. 

We show that the color-evaporation model predicts an inclusive branching ratio of Z 
into prompt ■i/' in agreement with the LEP results [|,^. The dominant process is Z — > ccqq. 
It predicts an energy spectrum which is very different from the one predicted by the color- 
singlet model as well as a substantially larger branching ratio of Z into charmonium. 

II. THE MODEL AND ITS PARAMETERS 

The color-evaporation approach, which actually predates the color-singlet approach, 
quantitatively describes all charmonium photo- and hadropro duct ion data |lTO| , p!T| . The 
model simply states that charmonium production is described by the same dynamics as Dt) 
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production, i.e., by the formation of a colored cc pair. Rather than imposing that the cc 
pair is in a color-singlet state in the short- distance perturbative diagrams, it is argued that 
the appearance of color-singlet asymptotic states solely depends on the outcome of large- 
distance fluctuations of quarks and gluons. These large-distance fluctuations are probably 
complex enough for the occupation of different color states to approximately respect sta- 
tistical counting. In other words, the formation of color-singlet states is a nonperturbative 
phenomenon. In fact, it does not seem logical to enforce the color-singlet property of the cc 
pair at short distances, given that there is an infinite time for soft gluons to readjust the color 
of the pair before it appears as an asymptotic ip, Xc or, alternatively, DD state. It is indeed 
hard to imagine that a color-singlet state formed at a range m^^ automatically survives to 
form a ip. This formalism was proposed almost twenty years ago and subsequently 
abandoned for no good reason. 

In the color-evaporation model the sum of the cross sections of all onium and open charm 
states is described by 

0-onium = 7: -7-^ , (1) 

9 J2mc dm 

and 

8 dace f , dace /„x 
c^open = 7: / dm — h / dm — — , (2) 

9 J2mc dm J2mD dm 

where the cross section for producing heavy quarks, acc, is computed perturbatively, irre- 
spective of the color of the cc pair, and m is the invariant mass of the cc pair. The coefficients 
I and I represent the statistical probabilities that the 3x3 charm pair is asymptotically in 



a singlet or octet state 10 



The color-evaporation model assumes a factorization of the production of the cc pair, 
which is perturbative and process dependent, and the materialization of this pair into a 
charmonium state by a mechanism that is nonperturbative and process independent. This 
assumption is reasonable given that the characteristic time scales of the two processes are 
very different: the time scale for the production of the pair is the inverse of the heavy- 
quark mass, while the formation of the bound state is associated to the time scale 1/Aqcd- 



Comparison with the data requires knowledge of the fraction of produced onium states 
that materiahze as ip^s, i.e., 

(^ip = Pip'^onium ; (3) 

where is assumed to be a constant. This assumption is in agreement with the low-energy 



data Notice that a single nonperturbative factor p describes a given charmonium 

state, regardless of the spin and orbital angular momentum of the charm pair. Therefore 
the color-evaporation model is more economical in nonperturbative parameters than the 
color-octet mechanism. 



In Ref. |]rT|, we determined the factor p^ from an analysis of photoproduction of char- 
monium. Using this value we were able to accommodate all data on the hadroproduction 
of charmonium without introducing any new parameters. This is a very remarkable result 
given that the subprocess responsible for the charmonium hadroproduction changes from qq 
fusion to gg fusion as the center-of-mass energy is increased. The fragmentation factor p^ 
turned out to be 0.50 (0.43) when we adopted the GRV 94 HO (MRS-A) parton distribution 
functions for the proton ( JT^). 



III. PROMPT DECAY OF Z INTO CHARMONIUM 

In the color-evaporation model the width for inclusive Z decay into prompt charmonium 

is: 

T{Z prompt charmonium) = g y dm > (4) 

where f is the partonic width for producing a cc pair. In order to obtain the partial width 
into a specific charmonium state we multiply the above expression by the appropriate frag- 



mentation factor p, which was determined from charmonium photoproduction data pO| , p 
Notice that the predictions for the Z decay into charmonium are parameter- free. Therefore 
the production of charmonium at the Z pole is a clean test of the validity of the color- 
evaporation model. 



We have evaluated all the tree- level partonic amplitudes using the package MAD GRAPH 
||16|| . The leading-order process in ag is Z ^ ccg, which leads to the production of a 
charmonium state and a hard jet. The partial width Z ^ tpg is only 6 x 10^^ GeV, 
for nic = 1.45 GeV, = 0.50, and as{2mc) = 0.235. It is small because the virtual quark 
propagator suppresses the amplitude by a factor of the order rric/mz- This situation persists 
even when taking into account the next-to-leading order corrections. Therefore, this decay 
mode is too small to be seen at LEP I. 

Although formally of higher order in a^, the dominant process for the inclusive decay 
of the Z into charmonium is Z ^ ccqq, where q = u, d, s, c, and b. We show in Fig. [l| 
the Feynman diagrams contributing to this process. Table I contains the color-evaporation- 
model predictions for the Z partial width into the different partonic final states [jl^. In 
order to access the uncertainties in the determination of the fragmentation factor from 
the photo-production of ip data, we considered the values p^ = 0.50 and 0.43 which were 
obtained from ip photoproduction using the GRV 94 HO structure functions with rric = 1.45 
GeV and MRS-A distribution functions with = 1.43 GeV, respectively [|n|]. Notice that 
the branching fraction of Z into prompt ip is (1.7-1.8) x 10~^. This is to be contrasted with 
the color-singlet model which predicts a branching fraction for direct ip in Z decay of the 



order 3 x 10"^ ||T8|. The color-evaporation model leads to a branching fraction larger by 
almost an order of magnitude. 

Using vertex detectors the LEP collaborations have recently been able to separate the 
prompt production of ip from those originating from 6-hadron decays The OPAL 

collaboration reports that 

B{Z prompt i/j + X) = (1.9 ± 0.7 ± 0.5 ± 0.5) x 10"^ . 

Therefore, the color-evaporation-model prediction is in excellent agreement with the LEP 
results, illustrating how this approach gives a complete picture of the charmonium produc- 
tion in hadron-hadron, 7-hadron, and Z decays. At this point we should point out that the 
color-octet model predicts a branching ratio Z ^ ip + X to he 2.9 x 10^^ which is also 



in agreement with the experimental data. 

In Fig. we show the energy distribution [dT/dz with z = 2E.^/Mz) of prompt if) 
for the dominant process Z — > ccgg, where we summed over all quark flavors. The band 
between the dashed curves is the distribution predicted by the color-evaporation approach, 
when we consider the different values for given above in order to access the theoretical 
uncertainties. For comparison, we show the same distribution for the dominant process 
[Z i'qq) computed in the context of the color-octet (singlet) model represented by the 
upper (lower) solid curve |ll9 |. 

As we can see from Fig. H, the color-evaporation distribution is softer than the one 
predicted by the color-singlet model; it dominates however over the entire allowed range 
for z. Moreover, the color-evaporation mechanism and the color-octet model predict very 
similar distributions. Taking into account the theoretical uncertainties in the determination 
of the nonperturbative parameters describing these models, it is impossible to distinguish 
between them on the basis of this distribution. Therefore the measurement of the energy 
spectrum of the prompt ip is not a distinctive signature of the color-octet model, contrary 
to what is claimed in Ref. |T^. In the OPAL analysis of the production of prompt ip^s [Q, 
it is shown that the color-octet model describes the momentum distribution of prompt ip. 
Since this model and the color-evaporation mechanism have similar predictions for this 
distribution, we can conclude that this last model also describes the available data. 



IV. CONCLUSIONS 

The predictions of the color-evaporation model for the production of prompt ip on Z 
decays have no free parameters once we use the photoproduction of charmonium to con- 
strain the nonperturbative parameters of the model. We showed that this approach gives 
a good description of the available data on the production of prompt at LEP I. Taking 
into account the success of this model to describe the photo- and hadropro duct ion of char- 
monium, we can conclude that this model gives a robust and simple parametrization of all 
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charmonium physics. 

We would like to stress that the color-evaporation model has the same degree of success as 
the color-octet mechanism; however, the number of free parameters in the color-evaporation 
approach is smaller. The ip^s produced through the color-evaporation mechanism are basi- 
cally unpolarized since the polarization information is lost because of the multiple soft-gluon 
exchanges |^0|. On the other hand, the (non)polarization of ip is hard to explain in the frame- 



work of the color-octet model plip2|. Therefore, the measurement of the polarization of the 
produced charmonium may very well be a tool to discriminate between these competing 
descriptions. 
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TABLES 

TABLE I. Partial decay widths for Z ^ tp + X. The results were obtained using = 0.50 
[0.43], rric = 1.45 [1.43] GeV, and as(2mc) = 0.235 [0.236]. 



X 


r(Z ^ V + A) (GeV) 


Br((Z^V + A) 


9 


6.0 10-'^ [5.6 10-'^] 


2.4 10-^ [2.3 10-'^] 


uu, dd, ss 


2.8 10-4 [2.6 10-4] 


1.1 10-4 [1.0 10-4] 


cc 


7.8 10-5 [7.3 10-5] 


3.1 10-5 [2.9 10-5] 


bb 


8.1 10-' [7.9 10-''] 


3.4 10-5 [3.2 10-5] 


qq (all q) 


4.4 10-4 [4.1 10-4] 


1.8 10-4 [1.7 10-4] 
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FIGURES 
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FIG. 2. Energy spectrum dT/dz of the V' from the dominant partonic process Z — > ccgg. The 
upper (lower) solid curve stands for the prediction of the color octet (singlet) model according to 
Ref. II^. The upper (lower) dashed curves is the distribution predicted by the color evaporation 
approach when we take = 0.50 (0.43), nic = 1.45 (1.43) GeV, and as{2mc) = 0.235 (0.236). 
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